Diffusion-free ultrafast carrier dynamics in silicon nano-pillars 
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We have investigated ultrafast carriers dynamics in crystalline silicon nano-pillars structure using 
a pump-probe reflectivity method with 800 nm, 150 fs laser pulses and fluence in the range of 
17 -j- 170 mj/cm 2 . Dimensions of the structure allow us to eliminate contribution from the diffusion 
process to the relaxation dynamics of the excited carriers. Strong intensity-dependent time-resolved 
reflectivity change, AR(rd), was monitored in the submelting regime. At strong pumping AR(rd) 
can reach up to 8%, a several times higher than for a bulk silicon. From the measurements we 
deduced recombination time of up to 10 ns at carriers density of ~ 1 x 10 20 cm" 3 , while electron- 
phonon interaction occurs during 350-400 fs and it is independent of the concentration. 

PACS numbers: 71.35.Ee, 72.20.Jv, 73.63.-b, 78.47.J-, 78.67.-n 



Since the femtosecond laser technology became avail- 
able to a wide range of research areas, the subpicosecond 
carriers dynamics in semiconductors has been subjected 
to intensive theoretical [l[ and experimental 0, 01, HI re- 
search in bulk materials. In recent years considerable 
attention has been devoted to the development of appli- 
cations based on Si nanostructures (NS). For example, 
Si field emitters based on nano- wires numerous ap- 
plications of NS in photonics: electroluminescence diode, 
microcavity and waveguide; electronics: single-electron 
device and nonvolatile memory Q . An understanding of 
carriers dynamics is crucial for the research and develop- 
ment of such devices. 

Transient measurements of reflectivity and transmis- 
sivity, provided by the pump-probe technique, allow tem- 
poral monitoring of changes in the optical properties of 
materials contributed by excited electrons (and holes) dy- 
namics. The main parameters governing carriers dynam- 
ics are: carrier-carrier scattering, occurring on a few fem- 
toseconds time scale Q; carrier-phonon thermalization 
lasting a few hundred femtoseconds § ; and electron-hole 
(e-h) surface [9( and bulk [lfj recombinations typically 
characterized by much longer time scales [TJ El ■ 

The drawback estimating of these parameters from the 
optical measurements is a carrier diffusion from the ex- 
cited area or 'probing spot', especially for high density 
e-h plasma, N e _ h > 10 19 cm" 3 p], El- This compli- 
cates interpretation of results owing to a convolution be- 
tween diffusion and various carrier-carrier and carrier- 
phonon interactions With the development of nan- 
otechnology it has become possible to construct a struc- 
ture where spatial confinement suppresses the diffusion 
process though optical properties remain similar to those 
of the bulk. Such a structure has to be smaller than the 
diffusion length (estimated as V^t^Djv ~ 1 — 100 /im, with 
the diffusion coefficient Dn — 10 cm 2 /s [111] ) but larger 
than electron wavelength to avoid effects of the quantum 



confinement. Thus structures with dimensions of about 
100 nm are suitable for the purpose. 

We present ultrafast transient reflectivity measure- 
ments of a crystalline Si nano-pillars (NP) structure 
where individual Si NP are sufficiently separated from 
the bulk substrate and each other. Thus, pump-probe 
measurements of an ensemble of NP represent individual 
properties of a single pillar free of contributions from its 
surrounding. In this Letter, we show: (i) carrier dynam- 
ics below and at the melting threshold (ii) dependence 
of the decay constants on the carriers concentration (iii) 
linear response to the excitation and consequently (iv) 
utilization of the Drude model. We also observed that 
the absolute change in the reflectivity is up to one order 
of magnitude larger than similar changes observed in bulk 
Si, and the carriers recombination time is concentration- 
dependent occurring on a timescale of a few nanoseconds 
while thermalization time is about a few hundred fem- 
toseconds independent of the carriers concentration. 




* Corresponding author: dimitri@nprl.ph.bham.ac.uk 



FIG. 1: SEM pictures of silicon nanopillars made by imprint 
method on Silicon-on-Insulator (SOI) wafer. 50° tilted view 
with resolution of 1 fim. Inset: 45° tilted view with resolution 
of 200 nm, pillar size estimation: 70 x 200 nm 
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The femtosecond pump-probe reflectivity measure- 
ments were performed by use of Thsapphire laser oscilla- 
tor amplified up to 1 mJ per pulse at 800 nm and rep rate 
of 1 kHz. Second harmonic generation of BBO crystal 
was used to determine the pulse duration of t p = 150 fs of 
nearly Gaussian shape. The pump and probe beams were 
p- and s-polarized respectively. Both beams had been 
nearly normal incident on to the sample (with a small 
skew angle with respect to each other) focused to a round 
spot with a diameter of 70 /im. Careful overlap of the 
pump and probe was confirmed prior to each measure- 
ment by a CCD camera equipped with 100 x magnifying 
objective. Probe ffuence was 200 /iJ/cm 2 while the pump 
fluence varied in the region of 17 4- 170 mJ/cm 2 . The re- 
flected probe pulse passed through linear s-polarizer, in 
order to reduce scattered p-polarized pump background, 
and was detected by Si detector. Another Si detector was 
used to monitor pulses intensity reflected from a glass 
wedge located before the sample. Signals of both detec- 
tors were averaged by a two-channel digital oscilloscope, 
the signal reflected from the sample was normalised to 
the monitored laser intensity signal. The remaining ex- 
perimental details can be found elsewhere [lij ■ Using this 
technique we are able to resolve the reflectivity changes 
of AR/Rq « 0.01 (defined as AR = R(r d ) - R , where 
Rq is the absolute reflectivity and Td is the delay time 
between the pump and probe). 

A nano-imprint assistant dry-etching procedure has 
been employed to fabricate Si NP on Silicon-on-Insulator 
(SOI) structure by Lamontagne et al [l5[. The pillars 
are separated by ~ 125 nm thick S1O2 layer from the 
bulk Si substrate. The density of the NP in the nano- 
mask is 2.5 x 10 9 cm -2 . Fig. [T] shows Scanning Electron 
Microscope (SEM) pictures of the sample with various 
resolutions. The highly-ordered structure of the NP is 
well resolved. The following average dimensions of the 
NP have been observed: diameter of ~ 70 nm and height 
of ~ 200 nm. The inset in Fig. [3] shows the absolute re- 
flection, Rq, measurement in the vicinity of A = 800 nm. 
According to this result, R = 0.075 ± 0.01 at A = 800 
nm, roughly four times lower than Rq of a crystalline Si. 
Similar value, of decreased Rq, was observed by G. R. 
Lin et al. [16| at an even broader spectrum on similar 
structures. 

In the scope of the Drude model with neglected damp- 
ing, the free carriers plasma refractive index is n p = 

noyl — w 2 /cjg. Here uq and n p are the refractive in- 
dices of the bulk at the ground state and modified by 
the contribution from e-h plasma respectively, loq is the 
probing frequency while the plasma frequency denoted 
as uj p — N e -he 2 / €om e f f ■ N e -h represents e-h plasma 
density, m e ff — {l/m e + l/mh)^ 1 is the reduced effective 
mass of the e-h pair, eo and e carry their usual meanings. 
Neglecting non-linear response, the relation between the 
change in the reflectivity and in the refractive index 
(An = n p — no) at normal incidence can be expressed 
as: AR/Rq = AAn/{n 2 - 1). Thus, AR cx An cx N e -h 



allowing extraction of N e -h from the AR/Rq experimen- 
tal data. 
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FIG. 2: (Color online) Transient reflectivity of silicon nanopil- 
lars as a function of the delay time at various fluences values 
indicated in F m units. Upper panel: in the vicinity of melt- 
ing fluence, F m - Lower panel: below melting fluence. Some 
curves are shown on both panels for easier comparison. 
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-linear fit, L = 1 ±0.05 
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FIG. 3: (Color online) Normalized change of reflectivity ver- 
sus fluence at delay time Td — 3 ps. Linear fit represents 
an order of the photoexcitation, L — 1 ± 0.05. Note: supra- 
linear scale. Inset: absolute reflection, Ro in the vicinity of 
A = 800 nm. 

The time resolved AR/Rq versus the delay time, Td, 
recorded on the Si NP at different fluences in the range 
of 0.1441 F m (where F m = 170 mJ/cm 2 , is the threshold 
melting fluence of Si surface [|[) is shown in Fig. O The 
upper panel shows that at the melting fluence, F — 1 F m , 
AR/Rq is predominantly positive. The initial reflectivity 
increase is attributed to the metal-like optical response 
contributed mainly by high-density e-h plasma which de- 
termines n p to be purely imaginary. After ~1 ps R[rd) 
decreases owing to carriers' collision-enhanced e-h re- 
combination process. The rest of the curves, recorded 
at the fluences below the melting point, present nega- 
tive AR/Rq. Below the melting fluence a minimum of 
AR/ Rq is reached after delay time of 1-2 ps. The largest 
change in the reflectivity, AR/Rq = 8%, is observed at 
F = 0.94 F m . For comparison, similar measurements 
near the threshold fluence on a crystalline Si wafer pro- 
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vided negative AR/Ro of ~ 1.5% (not shown). Fig. [3] 
shows a supralinear plot of \AR/Rq\ versus F at r d = 3 
ps. From the plot, the coefficient of non-linearity, L, has 
been estimated as 1±0.05, demonstrating primary linear 
response of the pillars to the excitation. 
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FIG. 4: (Color online) The extracted decay times vs plasma 
density, N e -h- Blue curve: the electron-phonon relaxation 
time, T e - P h- Red curve: the electron-hole recombination time, 
tv. Inset, upper part: experimental AR/Ro vs Ta at F = 
0.6 Fm and free carriers model fit. Inset, lower panel: N e _h 
vs Td is extracted from AR/Ro in the scope of the free carriers 
absorption model. 

In order to analyze the fast and the slow elec- 
tron dynamics, we have adopted the two-exponential 
model presented by Sabbah et al. [ItJ (see Eq's 6-8), 
here we note that for fitting we used AR(Td)/Ro = 
AR(0+)/R o x [erf(r d /r p ) + l]/2 x [1 + {mf™ K /m™f - 
l)exp(-r d /T e _ ph )]- 1 x [1 + Cexp(-r d /r r )]/(l + C), 
where C is a fitting parameter. According to the model 
the fast component describes electron-phonon interaction 
and the slow represents e-h recombination. On the exper- 
imental curves, below the melting fluence, the fast com- 



ponent is responsible for the initial decrease in the reflec- 
tivity, R(t,i), while the slow reflects its long recovery. The 
upper inset of Fig. 0] shows representative fitting of the 
model to the experimental AR/Ro at F = 0.6 F m . This 
experimental data was taken with higher resolution than 
that presented in Fig. O Similar analysis has been done 
for the decay curves measured at different fluences and 
the decay constants have been extracted. Fig. [5] presents 
these constants as a function of the initial N e _h which has 
been estimated assuming linear absorption with the co- 
efficient a = 1000 cm -1 and i?o is taken from the inset of 
Fig. [31 N e -h = a(l — R )F/E p h, where E p h is the photon 
excitation energy. The long decay time, r r , describes the 
e-h recombination time. It is found to depend strongly on 
the plasma density, decreasing from 10 ns down to 300 ps, 
at the density range of N e _ h = (0.8- 6) x 10 20 cm" 3 . The 
lower inset of Fig. [4] shows stability of the carriers pop- 
ulation (after initial change) on the timescale of tens of 
picoseconds at a fluence of F = 0.6 F m (N e _h = 3.6 x 10 19 
cm -3 ). The short decay, T e _ p /,, is found to be indepen- 
dent of N e _h and estimated to be ~ 350 fs. 

In conclusion, we have studied carriers dynamics in 
the excited crystalline Si NP. It is clearly shown that 
under strong diffusion suppression and at near submelt- 
ing regime, the recombination time of the e-h plasma is 
much longer than reported for a bulk material: by UV 
excitation [ll[ and diffusion dominated behavior 0] . The 
thermalization time, however, remains unaffected. The 
initial drop of the reflectivity is up to 8% which is a few 
times higher than for the bulk silicon. 

The authors appreciate the help of D C Sirica Ltd with 
the samples provided and thank E. Baskin for fruitful 
discussions during this work. We thank M. Y. Song for 
the SEM analysis. We acknowledge the financial support 
of AWM "Science City" network and UK EPSRC. 



[1] B. E. Sernelius, Phys. Rev. B 40; S. Jallepalli, M. 

Rashed, W. K. Shih, C. M. Maziar, and A. F. Tasch, 

J. Appl. Phys. 81, 2250 (1997). 
[2] C. V. Shank, R. Yen and C. Hirlimann, Phys. Rev. Lett. 

50, 454 (1983). 
[3] K. Sokolowski-Tinten and D. von der Linde, Phys. Rev. 

B 61, 2643 (2000). 
[4] A. Othonos J. Appl. Phys. 83, 1789 (1998). 
[5] Frederick C. K. Au, K. W. Wong, Y. H. Tang, Y. F. 

Zhang, I. Bello, and S. T. Lee, Appl. Phys. Lett. 75, 

1700 (1999). 

[6] N. Koshida, "Device Applications of Silicon Nanocrystals 
and Nanostructures" , Series: Nanostructure Science and 
Technology, Springer, 2009. 

[7] R. Buhleier, G. Lupke, G. Marowsky, Z. Gogolak, and J. 
Kuhl, Phys. Rev. B 50, 2425 (1994). 

[8] F. E. Doany and D. E. Grischkowsky, Appl. Phys. Lett. 
52, 36 (1988); W. Kutt, A. Esser, K. Seibert, U. Lemmer, 
and H. Kurz, Proc. SPIE 1268, 154 (1990). 

[9] A. J. Sabbah and D. M. Riffe, J. Appl. Phys. 88, 6954 



(2000). 

[10] A. Esser, W. Kutt, M. Strahnen, G. Maidhorn, and H. 
Kurz, Appl. Surf. Sci. 46, 446 (1990); F. E. Doany, D. E. 
Grischkowsky, and C. C. Chi, Appl. Phys. Lett. 50, 460 
(1987). 

[11] T. Tanaka, A. Harata and T. Sawada, Appl. Phys. Lett. 

82, 4033 (1997). 
[12] J. P. Woerdman Philips Res. Rep. 7, 12 (1971). 
[13] C. M. Li, T. Sjodin, and H. L. Dai, Phys. Rev. B 56, 

15252 (1997). 

[14] A. Kaplan, M. Lenner, and R. E. Palmer, Phys. Rev. B 

76, 073401 (2007). 
[15] B. Lamontagne et al, Department of Engineering Physics, 

Ecole Polytechnique, Montral, QC H3T 1J4, Canada. 
[16] G. R. Lin, Y. C. Chang, E. S. Liu, H. C. Kuo and H. S. 

Lin, Appl. Phys. Lett. 90, 1819231 (2007). 
[17] A. J. Sabbah and D. M. Riffe, Phys. Rev. B 66, 165217 

(2002). 



